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The Chinese Academy of Science has launched a thorium-based molten-salt reactor (TMSR) research project 
with a mission to research and develop a fission energy system of the fourth generation. The TMSR project 
intends to construct a liquid fuel molten-salt reactor (TMSR-LF), which uses fluoride salt as both the fuel and 
coolant, and a solid fuel molten-salt reactor (TMSR-SF), which uses fluoride salt as coolant and TRISO fuel. 
An optimized 2 MWth TMSR-LF has been designed to solve major technological challenges in the Th-U fuel 
cycle. Preliminary conceptual shielding design has also been performed to develop bulk shielding. In this 
study, the radiation dose and temperature distribution of the shielding bulk due to the core were simulated 
and analyzed by performing Monte Carlo simulations and computational fluid dynamics (CFD) analysis. The 
MCNP calculated dose rate and neutron and gamma spectra indicate that the total dose rate due to the core at 
the external surface of the concrete wall was 1.91 uSv/h in the radial direction, 1.16 uSv/h above and 1.33 uSv/h 
below the bulk shielding. All the radiation dose rates due to the core were below the design criteria. Thermal 
analysis results show that the temperature at the outermost surface of the bulk shielding was 333.86 K, which 
was below the required limit value. The results indicate that the designed bulk shielding satisfies the radiation 
shielding requirements for the 2 MWth TMSR-LF. 
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I. INTRODUCTION 


In 2011, the Chinese Academy of Sciences (CAS) 
launched a thorium-based molten-salt reactor (TMSR) relat- 
ed research project with a mission to research and develop a 
new fission energy system by 2020 [1]. The Molten Salt Re- 
actor (MSR) [2] is a class of nuclear fission reactors using 
molten salt as the fission and fissile material, as well as the 
coolant [3]. Unlike solid-fuel reactors, the fuel is dissolved 
in a fluoride salt and pumped into the reactor core and the 
heat is absorbed and removed by the molten salt itself in M- 
SRs [4]. MSR was selected as one of the candidates for the 
future development of nuclear energy by the Generation IV 
Forum (GIF) in 2002 [5] for its diverse advantages, such as 
good neutron economy, inherent safety, online reprocessing, 
less radioactive waste and nuclear nonproliferation [6, 7]. 

Since 2005, the GIF has focused on MSR from thermal 
spectrum molten salt reactor. Many works have been conduct- 
ed to investigate the radiation safety of the MSR. Elsheikh [8] 
presented the features of different safety characteristics of M- 
SR power plant and assessment in comparison to other sol- 
id fueled LWRs. Merk and Konheiser [9] applied the un- 
structured mesh neutron transport code HELIOS to determine 
the real neutron flux and power distribution for an advanced 
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molten salt fast reactor design. Zhu et al. [10] develope- 
d a one-dimensional neutron diffusion and transient analysis 
code to simulate neutron diffusion and transient in an MSR. 
In all the studies, the effect of the operation temperature was 
not taken into account. The TMSR project intends to con- 
struct two types of MSRs: one is a liquid fuel molten-salt 
reactor (TMSR-LF) which uses fluoride salt as both the fuel 
and coolant; the other one is a solid fuel molten-salt reactor 
(TMSR-SF), which uses fluoride salt as a coolant and TRISO 
fuel. A preliminary conceptual design has been carried out 
to develop a 2 MWth liquid fuel molten-salt reactor. Primary 
shielding calculations and design have also been performed 
for the 2MWth TMSR-LF. The primary goal of radiation 
shielding is to protect personnel from radiation hazards. The 
minimum thickness of materials to decrease the dose rate to 
an acceptable level should be determined. Another key point 
is the high operation temperature of the molten salt fuel. The 
fuel was designed to enter the core of the 2 MWth TMSR- 
LF at 908.15 K. The high temperature around the core was 
not suitable for shielding materials. So a water-cooled, cool- 
ing fins-filled thermal shield was designed to remove heat and 
bring down the temperature. The performance of the designed 
bulk shielding at attenuating the dose rate and bringing down 
the temperature should be assessed to confirm the radiation 
shielding safety for the 2 MWth TMSR-LE. It is important to 
note that only the core of the 2 MWth TMSR-LF was consid- 
ered as the radiation source, the auxiliary components such 
as the pipes of the primary system, the heat exchanger (IHX), 
and the pump were not taken into account in the analysis. 
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Fig. 1. Configuration of the core and the bulk shielding of the 2 MWth TMSR-LF. 


TABLE 1. Parameters about the fuel and the structures of the designed bulk shielding 


General Structures ID (cm) OD (cm) Height (cm) Materials 
Fuel BeF,-LiF-ZrF,-UF, Core - 140 160 - 
compositions 65-29.2-5-0.8 mole% 
Uranium fractions si aa eaga Core can 140 146 166 Hastelloy-N 

0.3-35-0.3-64.4 atom% Fuel 146 156 206 
Density 2270 kg/m? Reflector 156 166 256 Graphite 
Mass flow rate 171.86 kg/s Vessel 166 192 262 Hastelloy-N 
Fuel inlet temperature 908.15 K Thermal shield 240 320 540 Stainless steel & water 
Fuel outlet temperature 936 K Reactor cell 730 740 670 Stainless steel 
containment 

Fuel viscosity p[Pas] = 0.1542¢3624/7 Shield tank 740 900 835 Magnetite concrete 
Fuel heat capacity 1976 J/(kg K) Concrete wall 900 1060 995 Ordinary concrete 
Fuel thermal conductivity 1.44 W/(m K) 
Graphite heat capacity 1758 J/(m K) 


Graphite thermal d[w/(m + K)] = 3763 x Ty °7 


conductivity 


The As Low As Reasonably Achievable (ALARA) prin- 
ciple was used as guidance for radiation protection in the 
2MWth TMSR-LFE. Furthermore, the national standard for 
radiation source safety [11] and the associated guidance were 
utilized in the shielding design and analysis. Taking into ac- 
count these guidelines, shielding analyses were performed for 
the entire bulk shielding. Simulation of the radiation dose 
distribution and the neutron spectra due to the core were car- 
ried out for the bulk shielding of the 2 MWth TMSR-LF by 
using the Monte Carlo method. Thermal analysis was also 
carried out to present the temperature distribution of the de- 
signed bulk shielding for this high temperature reactor. 


I. METHOD 
A. Description of the 2 MWth TMSR-LF and bulk shielding 


The 2MWth TMSR-LF considered in this work is a 
2MWth graphite-moderated thermal reactor making use of 


solutions of uranium in fluoride carrier salts. A previous para- 
metric study was performed to develop this optimized MSR 
based on the Molten-Salt Reactor Experiment (MSRE) [12], 
which was originally developed at Oak Ridge National Labo- 
ratory (ORNL) in the 1960s. 

Shielding calculations were carried out and shielding de- 
sign was presented for the 2MWth TMSR-LF previously. 
The designed bulk shielding consists of the thermal shield, the 
shield tank, and the concrete wall, as shown in Fig. 1. Table 1 
gives information about the 2 MWth TMSR-LF and the bulk 
shielding. The reactor core, which was located in the center of 
the bulk shielding, was a 140 cm (diameter) by 160 cm (high) 
graphite matrix structure made up of 513 graphite stringer- 
s. The graphite stringers were mounted in a vertical close- 
packed array, as shown in Fig. 1. Half-channels were ma- 
chined in the four faces of each stringer to form flow passages 
in the assembly. In addition, the core of the optimized MSR 
was surrounded by a graphite lump, which not only serves as 
the moderator, but also provides a neutron reflector to gain 
good neutron economics. The optimized MSR was designed 
to be fueled with BeF,-LiF-ZrF,-UF, (65-29.2-5-0.8 mol%), 
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with atomic fractions of uranium of 0.3% 7°4U, 35% 735U, 
64.4% ?°8U and 0.3% *3°U. 

The thermal shield was a water-cooled, cooling fins-filled 
stainless steel container completely surrounding the reactor 
core. The thermal shield was designed to reduce the neutron 
dose rate by a factor of 104 and attenuate the gamma dose rate 
by a factor of about 10° in the horizontal plane. As mentioned 
above, the thermal shield was required to absorb the heat from 
the core significantly. It was designed so that the difference in 
temperature at the internal and external surface of the thermal 
shield was about 150 K. 

The fuel pump and heat exchanger were located in a 
370 cm diameter by 670 cm high reactor cell. This reactor cell 
was surrounded by a magnetite concrete shield tank, which in 
turn sits within an ordinary concrete wall [13]. The reactor 
cell was not allowed to be entered under any circumstances 
after the reactor was operated at power because of the high 
radiation level. The total dose rate criteria for the radiation 
shielding of the 2MWth TMSR-LF was 2.5 uSv/h at areas 
outside the concrete wall. 

Due to the principle of the temperature limit of the shield- 
ing materials, it was required by the ASME that the maximum 
temperature of the concrete should not exceed 366 K during 
operation [14]. 


B. Method of analysis 


The Monte Carlo code MCNP5 was employed as a cal- 
culation tool for numerical radiation transport studies. M- 
CNP is a general-purpose Monte Carlo N-Particle code that 
can be used for neutron, photon, electron, or coupled neu- 
tron/photon/electron transport, including the capability to cal- 
culate eigenvalues for critical systems. For the shielding cal- 
culation of the 2 MWth TMSR-LF, the fuel was considered to 
be static in the core. The total neutrons from the core and the 
fission gamma rays were considered as the radiation source in 
shielding analysis. The fission gamma rays were determined 
by source term calculation according to the operation plan. A 
MCNP model of the core and its bulk shielding was devel- 
oped following the geometry structure and dimensions given 
in Table 1 for the Monte Carlo simulations. The neutron and 
gamma spectra averaged over the core were calculated, then 
the spectra of each shield were simulated and converted to 
the dose rate by employing the flux-to-dose rate conversion 
factors recommended by ICRP Publication 74 [15]. The da- 
ta libraries used in this calculation were processed by NJOY, 
based on the most recent ENDF/B-VII nuclear data available. 

The computational fluid dynamics (CFD) code FLUENT 
was applied to the thermal analysis for the 2MWth TMSR- 
LF and its bulk shielding. 3D model of the core, and shields 
were set up by employing the 3D mechanical computer-aided 
design (CAD) program of SolidWorks. GAMBIT, a general- 
purpose preprocessor for CFD analysis, was applied to gen- 
erate meshes and define the boundary types of the model. 
The meshed model was then imported to FLUNET to sim- 
ulate the temperature field of the 2 MWth TMSR-LF and the 
bulk shielding. Materials, such as water and steel, were con- 
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tained in the FLUENT database, but other materials, such as 
the fuel, the graphite, and the Hastelloy-N [16], still need to 
be defined by the user. Boundary conditions, such as the flow 
rate of the fuel and the fuel temperatures at the inlet and outlet 
also need to be defined. The parameters for the CFD simula- 
tion are listed in Table 1. The heat source was defined based 
on the MCNP5 calculated power density of the core through 
user defined functions (UDF) [17-19]. 


Ill. RESULTS AND DISCUSSION 
A. Radiation dose distribution 


The neutron and gamma spectra, as well as the power den- 
sity, averaged over the core were analyzed first. The output 
of the MCNP5 was defined by a number of tally specification 
cards (Fn card). The Monte Carlo program provides seven 
basic neutron tally types and six standard photon tallies. For 
the measurements of the neutron and gamma spectra, the tal- 
ly for flux averaged over a cell (tally F4) was used. For the 
measurements of power density, the tally for fission energy 
deposition averaged over a cell (tally F7) was used. And the 
KCODE was employed to specify the source of the 2 MWth 
TMSR-LF. 15550 cycles with 10° particles per cycle were 
run for the simulation, and 50 cycles were skipped before the 
tally accumulation began. 
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Fig. 2. Neutron and gamma spectra averaged over the core of the 
2 MWth TMSR-LF. 


The neutron and gamma spectra averaged over the core 
were analyzed and shown in Fig. 2. Notice that the neutron- 
s are spread widely from 107° MeV to about 20 MeV. The 
peak in the neutron spectrum where most neutrons appear in 
the 2MWth TMSR-LF was between 0.01 eV and 1eV. The 
gamma rays have a relatively narrow bandwidth where the en- 
ergy is spread from 1 keV to about 20 MeV. There were two 
sharp peaks in the gamma spectrum around 0.1 MeV. 

Figure 3 indicates that the power distribution was symmet- 
ric in the axial direction of the core and the peak appears at 
the middle of the axis. The peak of the power did not appear 
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Fig. 3. Power density of the 2 MWth TMSR-LF. 


f Neutron dose rate} | ! i 
i- - ~ - Photon dose rate 
NNI Cae EEE 
~ AN & | T i 
< a] a | 5] 
A E 2) 2 
= a | A E 
2 a | | Bo 
S E | O] 
a | Reactor cell Sg] | 
i) a} ~j i 
n IB | x4 ! 
a] fee [Seu 
a l ae 
Rue ds h) Cees: 

0 60 120 180 240 300 360 420 480 540 


Radial Distance(cm) 


—~ 
g 
Sat 


Shield tank 
Shield tank 


Concrete wall 


Concrete wall 


Reactor cell 


Neutron dose rate 
i- - - - Photon dose rate 


0 100 200 300 400 500 600 700 800 900 1000 


Axial Distance(cm) 


Fig. 4. Radiation dose rate distributions of the 2 MWth TMSR-LF in the (a) radial and (b) axial direction. 


at the center of the radial direction because there were control 
rods at that position, as seen in Fig. 1. 


A transport calculation was carried out to simulate the ra- 
diation dose distribution under the condition of the maximum 
power in Fig. 3 to obtain a conservative estimate. For the mea- 
surements of the neutron and gamma dose rate, the flux-to- 
dose rate conversion factors written in the Dose Energy Card 
(DE) and Dose Function Card (DF) were used. Figure 4(a) 
presents the MCNP5S calculation results of the radiation dose 
rate at the radial direction along the center line of the core. 
Notice that the magnetite concrete thickness of 80cm was 
sufficient to reduce the neutron dose rate to 1 uSv/h. Howev- 
er, an additional 80 cm of ordinary concrete may be required 
to attenuate the photon dose rate. The total dose rates due 
to the core at the external surface of the concrete wall were 
1.91 uSv/h. Figure 4(b) shows the radiation dose rate at the 
axial direction as about 40 cm from the center line of the core, 
at which position the peak of the power appears. The total 
dose rate due to the core was 1.16 uSv/h above and 1.33 uSv/h 
below the bulk shielding. 


It also can be noticed from Fig. 4 that the neutron dose rate 
due to the core was at the 10° level at the internal surface and 
the 10? level at the external surface of the thermal shield. The 
photon dose rate due to the core was at the 104 level at the 
internal surface of the thermal shield and at the 10 level at the 
external surface. The attenuate factors of the neutron dose 
rate and photon dose rate meet the design objective. 

Figure 5 shows the neutron and gamma spectra outside of 
each shield at the radial direction along the center line of the 
core. The spectra indicate that the thermal shield has extraor- 
dinary attenuation in neutron flux and gamma flux. 


B. Temperature distribution 


As mentioned above, the molten salt was designed to 
be pumped into the core at 908.15 K, with a flow rate of 
171.86kg/s. Thermal analysis needs to be carried out to 
present the temperature distribution of the complex. The ther- 
mal analysis was displayed using the CFD code FLUENT. 
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The porous model [20, 21] was employed to simulate the core 
because the dimensions of the fuel passages were relatively 
small compared to the dimensions of the core. The power 
destiny distribution shown in Fig. 3 was loaded into the CFD 
model as the heat source through the UDF compiled in the C 
programming language. 

Figure 6 gives the temperature profile of the bulk shield- 
ing at the radial direction along the center line of the core, 
with the flow rate of the cooling water in the thermal shield at 
6.3 kg/s. Notice that the temperature difference between the 
internal surface and external surface of the thermal shield was 
about 150 K, which meets the design objective. The tempera- 
ture at the internal surface of the concrete wall was 333.86 K, 
which was below the required limit value. 


IV. CONCLUSION 


The radiation dose distribution and temperature distribu- 
tion of the 2 MWth TMSR-LF for the proposed TMSR project 
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was simulated and analyzed by performing Monte Carlo sim- 
ulations and CFD analysis. The neutron and gamma spectra 
from each shield were first calculated by MCNPS5S, then cal- 
culations were performed to present the total radiation dose 
distribution, based on these spectra, by using the flux-to-dose 
rate conversion factors. The calculated results indicate that 
the total dose rate due to the core at the external surface of 
the concrete wall was 1.91 uSv/h in the radial direction and 
1.16 uSv/h above and 1.33 uSv/h below the bulk shielding. 
The neutron and gamma spectra show that the thermal shield, 
composed of stainless steel fins and cooling water, attenuates 
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